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Abstract. The limits imposed on diffusive shock acceleration 
by upstream ion-neutral Alfven wave damping, and by ionisa- 
tion and Coulomb losses of low energy particles, are calculated. 
Analytic solutions are given for the steady upstream wave ex- 
citation problem with ion-neutral damping and the resulting 
escaping upstream flux calculated. The time dependent prob- 
lem is discussed and numerical solutions presented. Finally the 
significance of these results for possible observational tests of 
shock acceleration in supernova remnants is discussed. 
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1. Introduction 

Diffusive shock acceleration is generally believed to be an 
important astrophysical mechanism for pr oduci ng high en- 
ergy char ged p articles (Blandford & Eichler |l987|; Berezhko & 
Krymsky 



1988: Jones & Ellison 



1991) 



In particular, operat- 
ing at the strong shocks of Galactic supernova remnants, it 
is thought to be responsible for producing the Galactic cos- 
mic rays i.e., those of energy less than about 10 14 eV. However 
this hypothesis has not been convincingly confirmed by any 
direct observational test. Possibilities which have been sug- 
gested include the observation of high energy gamma-ray emis- 
sion from supernova remnants, or detection of shock prec ursor 



structures via B almer line observations (Aharonian et al. 1994 
I)' 



1991 



Raymond 

A key element in the effective operation of the acceleration 
process is the resonant excitation of scattering waves by the 
accelerated particle pressure gradient ahead of the accelerating 
shock resulting in much smaller values of the particle diffusion 
coefficient near the shock t han in the general medium. This 
was first emphasised by Bell (1978) who gave an analysis of the 
process in the steady state and also pointed out the importance 
of ion-neutral friction in damping the waves and quenching the 
acceleration of high energy particles. 

Another important aspect of diffusive shock acceleration 
is that it can only operate in conjunction with an injection 
mechanism, which must accelerate particles directly out of the 



'thermal' plasma up to a velocity of several times the thermal 
speed. However in dense media low energy injected particles 
are subject to Coulomb and ionisation energy losses. The ac- 
celeration process has to be fast enough at energies somewhat 
above 'thermal' to compete with these collisional processes, 
otherwise the shock will only be able to accelerate the few am- 
bient pre-existing high energy particles. 

These effects obviously have important implications for ob- 
servational diagnostics such as the gamma ray luminosity of 
supernova remnants and, in the case of non-radiative shocks, 
the detailed line profiles of the faint Balmer line emission. The 
gamma-ray emission will only be easily detected if the shock 
is propagating into a relatively dense medium witha num- 
ber density in excess of 0.1cm -3 (Aharonian et al. 1994) and 
the Balmer diagnostics require that the medium into which 
the s hock is propagating be only about 90% ionised (Ray- 
mond 



1991 



Rather generally one can anticipate that most 
observational tests for cosmic ray acceleration in supernova 
remnants will require some significant matter density to yield 
a detectable signal. However too high a density may quench 
the acceleration and destroy the very effect one is trying to 
detect. 

In this paper we first develop the theory of wave damping 
and excitation in cosmic ray shocks in more detail than has 
been done previously and derive corrected estimates for the 
energy at which this process quenches the acceleration. Then 
we consider the conditions under which collisional processes 
can suppress the natural shock injection process. Finally we 
discuss the implications of these results for the observability of 
cosmic ray acceleration effects in supernova remnants. 

2. Wave excitation and damping 

The model we study is identical to that considered by Bell 
although we use a slightly different notation. We ignore the 
reaction of the accelerated particles on the flow, and consider 
only their reaction on the scattering waves, which we take to 
be Alfven waves, moving at speed V along the magnetic field 
relative to a steady upstream flow of speed U and density p in 
the shock rest frame. Transmission of the waves through the 
shock and their damping downs tream has been considered by 
Achterberg and Blandford 



1986 
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and is not addressed here. For 
simplicity, we will consider the case of parallel shocks, in which 
the magnetic field lies along the normal to the shock surface. 
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Particles of momentum p, charge e and pitch fi (fj, = cos 9 
where 9 is the angle between the particle's momentum vec- 
tor and the mean magnetic field direction) resonantly interact 
with waves of spatial wavenumber along the field k — l/(p,r s ) 
where r g = p/(eB) is the gyroradius of the particle in the field 
of strength B (there is an implicit assumption here, that the 
waves are 'slow' and the particles 'fast'). Although there is a 
contribution to the scattering of particles of momentum p by 
waves on all spatial scales smaller than the particle gyroradius, 
it is clear that the bulk of the scattering is from waves close 
to r s in scale. In common with most discussions of this topic 
we assume that particles of momentum p, which from now on 
we take to mean particles in a (natural) logarithmic interval of 
momentum, interact only with waves in a logarithmic interval 
of k space where k and p are implicitly related by kp = O(eB). 
This is often called "sharpening the resonance." 

3. Equations of the model 

The accelerated particles diffuse against the flow so that the 
spatial and temporal evolution of the isotropic part of the par- 
ticle phase space density f{t,x,p) is given by the advection- 
diffusion equation 



dt ( ' dx dx 



with a diffusion coefficient k which we take, from quasi-linear 
theory, to have the form 



KB = 



(2) 



where T is the dimensionless resonant wave intensity per loga- 
rithmic interval of wavenumber and Kb is the so-called Bohm 
diffusion coefficient corresponding to a mean free path of order 
the particle gyroradius. The total wave energy density is given 
by 



( SB2 ) _ (BY 



2(i 



2p 



X(fc)dln(fc). 



(3) 



Physically all this says is that the larger the fluctuations SB 
imposed on the mean field B, the stronger the scattering, until 
when SB m B the field is totally disordered and particle tra- 
jectories lose coherence on scales as short as the gyroradius. 
Note that the advection is with the waves which we assume to 
be moving backwards at velocity V relative to the background 
plasma, and thus to be moving at velocity U — V relative to 
the shock (which we locate at x = 0). 

In scattering off the waves the particles do work at the rate 
VVP where VP is the accelerated particle pressure gradient. 
The resonance sharpening simplification allows us to apply this 
to a particular set of resonant particles and waves to give the 
wave energy equation, 



at 



+ (U- 



dx dx 



7 J 



(4) 



where I is the (dimensional) wave energy density, P the reso- 
nant particle pressure and 7 a damping coefficient. It is con- 
venient to write this in nondimensional form, and scale / to 
the background field energy density B 2 /2/io = piV 2 /2 and P 
to the ram pressure of the background flow of ionised plasma, 



piU , where pi is the ion mass density. We let V denote the 
dimensionless pressure per logarithmic momentum interval, 



Pit/ 2 3 jyP> 
so that the total accelerated particle pressure is 

T^//(P) dp = Pi [/ 2 / P{p)d\n{p) 



(•>) 



(6) 



Expressed in terms of the dimensionless quantities V and X we 
then obtain the fundamental pair of coupled equations 



dP lrr 
and 



v) dP = d_ f^B&P\ 
dx dx \ X dx J 



m di 2U 2 dp 
v) d-x = —^~^- 



(7) 



(8) 



We assume that the main process contributing to the wave 
damping is ion neutral friction. If the neutrals do not move 
coherently with the ions, the resulting wave damping rate 7 
has the value 



(1) 7 = n n {av) 



(9) 



where n n is the neutral number density and {av) is the mean of 
the collision velocity t imes the charge exchange cross section; 
Kulsrud & Cesarsky ( 1971) give the approximation 



{av} w 8.4 x 10" 



(— V 



(10) 



for temperatures T in the range 10 2 K to 10 5 K. The assump- 
tion that the neutrals do not participate in the coherent oscil- 
lations of the ions in the Alfven wave implies that the period 
of the wave must be short relative to the momentum transfer 
time scale from the ions to the neutrals. By assumption the 
wavelength is of order the gyroradius of a resonant particle, 
thus this condition implies 



< 



V rii {av} 



(11) 



where m is the ion number density. Numerically, this condition 
translates, for accelerated protons of energy E, to 



E 



lGeV 



< 8 1 



T 



10 4 K 



1 cm -3 



-3/2 



B 
l(j,G 



(12) 



At higher energies the resonant waves are of sufficiently long 
period for the neutrals to participate coherently in the wave 
motion. This will lead to decreased wave damping and a reduc- 
tion in the effective Alfven speed. The full dispersion relation 
for Alfven waves in a parti ally ionised medium is dis cusse d 
by Kulsrud & Pearce fll96S|) and also by Volk et al. fll98l|). 



However for our purposes what is required is not really the 
dispersion relation for freely propagating waves but the damp- 
ing of driven modes. An elementary calculation of the energy 
dissipation in a frictionally coupled system of ions and neutrals 
when the ions are forced at frequency u> gives 



OJ 2 + L0( 



U>i,n = ni, n {ov) 



(13) 
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Thus, as one would expect, at high frequencies the damping 
is constant and equal to ui n . However as u> — > u>i the damping 
rate drops to w n /2 and at lower frequencies decreases as uj 2 . 
Because the resonant wave frequency is inversely proportional 
to particle momentum p at sufficiently high particle energies 
the damping rate of the resonant waves will decrease as p~ 2 ■ 



4. Parameters 

Looking at this system we see that it explicitly contains only 
two dimensionless combinations of parameters. One is the ratio 
V/U which is the inverse Alfven Mach number of the shock. 
The other is jk^/U 2 which is essentially the ratio of the ac- 
celeration time scale to the damping time scale. However the 
boundary conditions on X and V introduce additional dimen- 
sionless parameters into the problem. 

The ratio V/U, for astrophysical shocks, is expected to lie 
in the range 10 -3 to 1. In fact the assumption made above, 
that the amplified waves are those propagating backwards rel- 
ative to the flow breaks down unless V < U. For most of the 
calculations reported here we have used V/U = 10 -2 . 

The dimensionless wave damping parameter due to ion- 
neutral friction has the approximate numerical value at ener- 
gies below the limit given by dl 2ft 



2^ = 25 2 ( U » } (-?—Y A (_M 
U 2 Vlcm-V U0 4 K7 VlTeVy 

-l 



cm 
B 
1 fj,G 



U 



10 3 kms- 



TeV, 

2 



(14) 



where n n is the neutral number density, T the temperature 
and E the accelerated particle energy (for relativistic protons). 
Because the acceleration time scale increases with energy the 
effects of wave damping become more significant as the energy 
increases to the value given by (jl^). At higher energies the 
damping rate decreases faster than the acceleration rate. Thus 
the energy given by ( p^ ) represents a critical threshold at which 
the wave damping effects are strongest. If the acceleration can 
reach this energy, then the ion-neutral wave damping alone 
cannot restrict further acceleration and the upper cut-off will 
be determined by other factors such as shock geometry or age. 

The dimensionless resonant particle pressure at the shock 
can be estimated by assuming that the shock acceleration is 
efficient, so that the total particle pressure is comparable to the 
ram pressure of the flow, and that the spectrum is roughly that 
predicted by test-particle theory for a strong shock, namely 
f(p) oc p~ 4 up to a maximum momentum p max (which may be 
of order 10 14 eV for supernova remnant shocks). This gives 



V 



ln(p max /rac) 



i 0.1-0.05 



(15) 



as the (steady state) value at x = 0, which is the right hand 
boundary condition on "P. In general one could have a popu- 
lation of particles advected into the shock from far upstream 
(a; — » — oo), corresponding to a finite far upstream value of V, 
however the physically interesting case is where all the parti- 
cles are produced at the shock so that the left hand boundary 
condition is V — *■ as x — > — oo (note that if I also tends to 
zero this is compatible with a nonvanishing escaping flux of 
particles upstream). 

For the resonant wave intensity, we have to specify a left 
hand boundary condition. One tempting possibility is to sup- 
pose that all the waves are created by the resonant particles, so 



that 1 — » as i - >—oo. However this is a rather singular limit 
as we will see in the next section. In general one has to sup- 
pose that some "seed" field of waves exists in the far upstream 
plasma and is advected into the shock where the backward 
propagating ones are then amplified by the resonant particles. 

5. Steady analytic solutions 

If we set the time derivates to zero in Eqs. (^) and (Eh, these 
can be linearised by introducing a new independent variable r 
similar to the "optical depth" variable used in radiative trans- 
fer: 



dx 



, U — V 

k(x') 

, s U - V 



(16) 



dx T(x 



where x~ is an arbitrary point in the flow. We then find 
dp d f, rr „,dP\ 



(U-V) 



dr dr 



((U-V) 



dr J 



and 



(U-V) 



dl _ 2lT_d£_ _ "/kb 

~d7 ~~ ~v~~d7 ~ u-v 



(17) 



(18) 



Equation (|17]) immediately yields a conserved quantity - the 
flux <f> of particles: 



{u -v)(v 



dV 



drj 



(19) 



and imposing the boundary condition V — V- at x = x- leads 
to the solution 



V 



(v- -+-) 
v u-v) 



+ [V 



u-v 1 v- u-v (20) 

This may be substituted into Eq. (|l^) , leading to the solution 

2U 2 [V-{U -V) -</>] , T 



2 — 1- 



KB 7 



(U-V) 



V(U - V) 2 



(e T -l). (21) 



Although it is in general both inconvenient and unneces- 
sary to evaluate V and T as functions of x, in the case of 
no damping 7 = 0, the procedure is straightforward. Defining 
a = 2U 2 [P-(U - V) - <P]/[{U - V) 2 V] and (3= (a-I-)(U- 
V) (x — x- )/kb , one finds 



a — T- 



so that 
V = 



+ 



aV(U -V)(a-I-) 



U-V 2U 2 {a-l-el 3 ) 



and 



1-e 



a-1- 



(22) 



(23) 



(24) 



a -X-eP 

These solutions diverge at a point x\ given by 

xi = x- + K B \n{a/l-)/[{a-l-.){U -V)] , (25) 
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which must therefore be chosen to lie in the downstream re- 
gion. In the special case a — 0, we recover the solution found 
by Lagage & Cesarsky (1983) by allowing x- — ► — oo, whilst 
keeping x\ finite: 



X = I_ < 1 — exp 



(U - V)l-(x - xi) 



which, in the limit X- 
solution found by Bell 



(26) 



0, reduces to the somewhat simpler 



1978) 



X = 



kb 



(U-V)(xi-x) 



(27) 



This is also a good approximation to the more general solution 
Eq. (|2^) for X ^> X- . We remark at this point that the trans- 
formation Eq. (jl6|) has eliminated the uninteresting solution 

V =constant, 1 = 0. Nevertheless, I_ — * is a singular limit 
of the full set (|) and (§). 

When damping is included, we note from Eq. (2lj) that 
there exists a point where the wave growt h is ju st balanced 
by damping and dX/dr = 0. Following Bell (1978), we assume 
that at this point the particles escape freely into the upstream 
plasma. As usual in the diffusion approximation, we can imple- 
ment this boundary condition by demanding that the density 

V vanish there. Let us choose x = X- and r = at the free es- 
cape boundary. Since the derivative of X there is by definition 
zero and we also demand V- = 0, the escaping particle flux is 
determined independently of X_ to be: 



kb7V 
' 2U 2 



(28) 



As pointed out by Bell, the escaping flux will steepen the spec- 
trum if it is comparable to the flux advected away downstream 
Thus, wave damping leads to a cut-off in the spectrum where 

7Kb V ^U-V 
U 2 2 



-Vo 



(29) 



wher e Vg is the value of V at the shock. Note h oweve r that 



Bell (1978), and following him Draine and McKee (1993), over- 
estimate the escaping flux by a factor c/(U — V) where c is 
the speed of light in vacuo. This error resulted from the use 
of the solution given in Eq. (j2^), in which the particle flux is 
identically zero everywhere, as a basis for the estimate of the 
escaping flux. 

6. Stability and time dependence 

If we examine the system of equations and think about the 
physics involved, it is clear that the wave excitation is strongest 
in those regions which already have an enhanced level of wave 
activity, and thus smaller values of the diffusion coefficient and 
steeper pressure gradients. It follows that the stability of the 
system is rather questionable; there is a clear possibility of 
some form of modulational instability. It is also possible that 
the steady solutions are not in fact the physically realistic ones. 
Before examining these questions numerically, it is interesting 
to look at the case of small high-frequency perturbations, which 
can be examined analytically. 

We begin by linearising about a smooth background solu- 
tion, V and X, on which is superimposed a small fluctuating 
component, V,X. The linearised equations are 



DV _ d_ ( KB &P _ ^Bf&P 

Dt ~ dx\ X dx X 2 dx 



(30) 



and 

DX _ 2U^_dP - 
~Dt ~ ~~dx ~ 7 

where 

D d , rr d 
D-t = 8t +{U ~ V) lTx 



(31) 



(32) 



denotes the advective derivative. We now assume that the fluc- 
tuations are of high spatial and temporal frequency and make 
a formal expansion in inverse powers of the wavenumber k, 



X = e ie{x ^ k~ n X n 



(33) 
(34) 



where 6(x,i) is the rapidly varying phase, 88 /dx = —k and 
d6/dt = uj = O(k). Substituting and collecting similar powers 
of k we obtain, to order k 2 and k, 



Po = 

V l = l^ET 

X x 

DX .2[/ 2 ~ - 



(35) 
(36) 

(37) 



together with the relation u = (U — V)k. Combining these 
results we get 



DXo _ 2U^_XodP _ ~ 
Dt V X dx 7 

which we can write in the form 



1 DX 
Jo Dt 



1 DX 

xUt' 



(38) 



(39) 



It follows that, to lowest order, small-scale perturbations in 
the resonant wave intensity are advected in at velocity U — V 
and grow (or decay) at exactly the same rate as the general 
background wave intensity. While this is not an instability, it 
does mean that the system "remembers" all the fluctuations in 
the initial seed wave field and is neutrally stable; in this lowest 
order linear high frequency analysis the relative amplitude of 
perturbations neither grows nor decays. Clearly a numerical 
study of non-linear finite size perturbations is desirable. 

There is no particular problem in developing a time- 
dependent numerical scheme for solving this system of equa- 
tions. However some thought has first to be given to a time- 
dependent right hand boundary condition for V . The problem 
arises because we are trying to look at one small part of the 
general process of shock acceleration in isolation. Essentially 
we have to model the rest of the shock acceleration process in 
a simple, but physically plausible way, in the boundary condi- 
tions we use. We assume that there is a steady flux, of particles 
into the momentum interval we are considering from accelera- 
tion at lower momenta. We are not interested in modelling the 
downstream scattering, so we simply assume that there is effi- 
cient scattering downstream and that the accelerated particles 
fill a downstream phase space "reservoir" of size 4ivp 3 L with 
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L ~ kb/W 2 where W 2 is the downstream advection velocity 
(we avoid discussion of whether this is the plasma velocity mi- 
nus the Alfven velocity). From this they can diffuse into the 
upstream region and they can be removed either by advection 
downstream, or by being accelerated to even higher energies. 

We need an equation for the resonant particle pressure, V, 
at the shock and in the downstream "reservoir" . The required 
equation is simply energy conservation as applied to the reser- 
voir. If we denote by Wi = U—V the advection velocity into the 
shock and by £ the energy density in resonant particles, then 
the energy flux in resonant particles into the reservoir from up- 
stream is Wi(£ + V) — n(d£/dx), the flux out downstream is 
Wi(£ + V), the flux in from lower momenta is £o(Wi — W 2 )/3 
where £0 is the energy density per logarithmic interval at lower 
energies and the flux out to higher energies is £(W\ — W 2 )/3. 
Thus 



T d£ 

1j — — 

dt 



Wl (£ + V)-^ 
2 ox 

-W 2 (£ + V)- 1 2 (£ - £o). 



(40) 



Note that, because we are only considering particles in a spe- 
cific momentum interval, this is not a two-fluid approximation. 
To simplify matters we assume that W 2 = Wi/4 and that we 
are dealing with relativistic particles where £ = 3V. The right 
hand boundary condition for V then becomes simply 



L ^ = WW > + 37>)-^ 
dt 2 Ox 



(41) 



where Vq is the value to which V tends if the solution ap- 
proaches the standard steady result (note that in a steady so- 
lution 



KB &P 

2 dx 



4W 2 V). 



(42) 



7. Numerical solutions 



It is not hard to solve the fundamental system of equations (J7J 
^ numerically. We have used a uniform grid and fixed time 
step with Ax/ At — (U — V) so that the advection is exact. 
The diffusion is treated fully implicitly giving a simple tridiag- 
onal system and the wave intensity is defined on a staggered 
mesh between each pair of V values. The wave damping is also 
treated implicitly so that the scheme is absolutely stable for 
all diffusion coefficients and damping rates. For more elabo- 
rate calculations there would be definite advantages in going 
to variable grid spacings and time steps, and perhaps a second 
order scheme (our scheme is only of first order). However for 
investigating stability there are distinct advantages in using a 
simple first order scheme with no advection errors. 

We first demonstrate that this program can reproduce the 
analytic results for steady solutions. In Figure [j] we show the 
solution attained at time 100 (in units of Kb/U 2 ) starting from 
V — and 2 = 0.1. The parameters have been chosen so that 
in Bell's solution ^2l\ ) the right-hand value of T is unity. The 
left-hand boundary conditions are zero net flux for V and fixed 
2 — 0.1 for the "seed" waves. In Fig ^| we show 1/2 over an 
extended range and for various upstream seed values. Note that 
with an upstream value of 0.01 the plot of 1/2 is an almost 
perfect straight line, as required by Bell's solution, and that 
the solutions for other seed values tend asymptotically to this 




Fig. 1. The numerical solution after time IOOke/l^ 2 for V/U = 0.01, 
7 = and "Po = 0.005. The dashed line is X and the solid line is 
V/Vo- Distance upstream is in units of kb/U. 




Fig. 2. The reciprocal of the resonant wave intensity, 1/1, after time 
200k b /U 2 in the region from —IOOk^/U to 0. The parameters are 
as in Fig. [t] except that the the left-hand value of X is successively 
set to 0.005, 0.01, 0.0125, 0.01667, 0.025 and 0.05. 



solution in the manner described by Lagage and Cesarsky's 
solution. 

In figure Fig. |^ we show the effect of including modest 
damping on the solutions of Fig. ^. As is clear physically, and 
obvious from the analysis of section 4, far upstream the damp- 
ing dominates. However the various solutions do converge to a 
single asymptotic solution near the shock. This is rather eas- 
ier to see in Fig. ^ where the final sections of the solutions of 
Fig. ^| are plotted in the style of Fig. ^. The left-hand bound- 
ary condition for these solutions was taken to be the constant 
escaping flux given by equation |28[ It is then easy to see that 
the right-hand boundary condition (|4l|), in the steady state, 
implies 



(43) 
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tion. Calculations with different periods and amplitudes yield 
similar results; although locally very different, the large-scale 
structure appears to be well represented by the steady solution 
corresponding to the "average" value of T. 



Fig. 3. The reciprocal of the resonant wave intensity, with 
modest damping after time 200kb/££? in the region from — 100kb/U 
to 0. The parameters are as in Fig. ti except that the damping rate 
7 = 0.1. 



V 1/ 



— BO -60 



which for the parameters used here gives V = 0.003 or V /Vo = 
0.6, as shown by the numerical solutions. We note that this ar- 
gument gives an approximate form for the spectrum at energies 
below the cut-off given by equation (E9) . 



Fig. 5. The reciprocal of the wave intensity with V/U = 0.01, 
7 = 0.1, Vq = 0.005 after time 200kb/C/ 2 with a sinusoidally mod- 
ulated and a steady left-hand boundary condition for I. The solid 
curve represents the solution with I = 0. 01(1 + 0. 9sint) as boundary 
condition, the dashed curve that with I = 0.01. 




-80 -60 -40 -20 



Fig. 4. The resonant wave intensity I (dashed lines) and V/Vo 
(solid lines) near the shock for the six solutions of Fig. til 

However it is most unlikely that the upstream "seed" field 
will be absolutely steady. In Fig. |B| we show the solution from 
Fig. ^ corresponding to a steady left-hand boundary condition 
of I = 0.01 and also that of the strongly modulated boundary 
condition, T — 0.01(1 + 0.9 sin t). More revealingly, in Fig. o, 
we plot the ratio of these two solutions. It is remarkable that, 
even at this very high level of modulation, the solution be- 
haves exactly as predicted by the analytic theory. Only very 
close to the shock is there some evidence of a slight increase 
in the perturbation amplitude. It is also noteworthy that the 
steady solution corresponding to the "average" value of T is 
a good approximation to the "average" time-dependent solu- 



Fig. 6. The ratio of the two solutions shown in Fig. Q 

In fact one has to be rather careful about exactly what 
sort of average is meant in the above statements. If one is spa- 
tially averaging along the particle gradient, the appropriate 
mean diffusion coefficient is the harmonic mean, correspond- 
ing to the simple arithmetic mean of the wave intensity. If one 
is averaging in the perpendicular plane, the appropriate mean 
is the arithmetic mean of the diffusion coefficient. More gen- 
erally one has to consider the geometry of the fluctuations as 
in percolation theory. It may be helpful to think of the electri- 
cal analogy of resistance networks connected in series and in 
parallel. To some extent, the fact that the arithmetic mean of 
the wave intensity generates a steady solution which is a very 
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good approximation to the "average" non-steady solution, is 
an artifact of the strictly one-dimensional model we are using. 

Although many extension or generalisations of diffusive 
shock acceleration have been considered, as far as we are 
aware there has been no work on shock acceleration with 
time- dependent scattering. If one considers the standard mi- 
croscopic arguments leading to the steady-state spectrum for 
test-particle acceleration, it is clear that the time-averaged 
spectrum at the shock will be the same power-law spectrum, 
independent of the details of the time-dependent scattering, 
as long as three conditions are met. All particles heading up- 
stream must return to the shock, which will be the case for 
negligible wave damping or if there exists sufficient upstream 
wave power. The angular distribution of the accelerated par- 
ticles at the shock front must be close to isotropic, which is 
generally required for the diffusive description to be applica- 
ble. And, crucially, the mean flux advected away downstream 
must equal the downstream advection velocity times the mean 
density at the shock. This last assumption appears open to 
question although it is in fact used in our approximate right- 
hand boundary condition. This may be an interesting way to 
produce non-standard spectral indices. 

8. Summary of wave damping effects 

In the preceding sections we have generalised the analytic the- 
ory of steady resonant wave excitation to include wave damp- 
ing and obtained an analytic formula for the resulting escaping 
particle flux. In reality the solutions will be time-dependent, 
but analytic arguments and numerical experiments show that 
"on average" the steady solutions can still give a good repre- 
sentation of the system (although this would probably change 
if more complicated non-linear processes were considered). 

The limit we derive for the maximum energy to which par- 
ticles can be accelerated before the escaping upstream flux 
kills the acceleratio n diffe rs substantially from that quoted in 
Draine and McKee (1993). Their formula (4.13) gives 
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where in their notation _E u g is the upper cut-off energy for 
a proton in units of GeV, v B 7 is the shock speed in units of 
10 7 cms _1 , no is the total number density and Xi the ionisa- 
tion fraction. Note that their dimensionless particle pressure 
parameter <f> is related to ours by 



= x% I Vodlnp 
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Our result ( |2S| ) says that if the upper cut-off energy E is 
determined by the ion-neutral wave damping in the upstream 
region, then numerically 
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where n n is the neutral number density and rii the ion density 
in the medium into which the shock is propagating. Our limit 
is a factor c/v B = 3x 10 3 /v B 7 times higher. In fact, as discussed 
above, the damping will become unimportant above energies 



given by (|l2|). The condition that the upper cut-off energy 
exceed that at which the neutrals start to coherently move 
with the ions is 
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If this condition is satisfied upstream ion-neutral wave damping 
places no restriction on shock acceleration. 

9. Collisional losses near injection 

At low energy, the most important loss processes for cosmic 
rays are ionisation and Coulomb losses. If the theory of diffusive 
acceleration at shocks is to be viable, it must overcome these 
for all energies above that of injection. The constraints that 
this imposes on the density, temperature and magnetic field 
of the ambient medium as well as on the energy of injection 
can be estimated simply by comparing the loss-rates with the 
acceleration rate. 

Standard treatments of diffusive acceleration give the ac- 
celeration rate in terms of the rate of change of momentum, 
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where the subscript refers to the upstream (1) or downstream 
(2) medium, k is the diffusion coefficient and u the fluid speed. 
However loss rates are normally discussed in terms of the par- 
ticle kinetic energy; defining the energy acceleration rate i^cc 
as the rate of change of kinetic energy divided by the kinetic 
energy we get 
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where r = ui/ita is the compression ratio of the shock, /3c is the 
speed of the particle and the Lorentz factor is 7 = (1 — /3 2 )~ 1 ^ 2 . 
The diffusion coefficient can be conveniently parameterised in 
terms of the Bohm value in a magnetic field B: 
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leading to 
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At non-relativistic energies the acceleration rate is proportional 
to /3~ 2 . 

Expressions for the l oss p rocesses have been given by 
Mannheim & Schlickeiser 1994. For Coulomb collisions with 
free electrons of temperature T = T4 x 10 4 K and density 
XiriH cm -3 , where x\ is the ionisation fraction and nn the total 
(ionised, neutral and molecular) density of hydrogen atoms, 
one finds 



ip 1 = 3.3 x 10" 16 XinH 
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where /3th c is the thermal speed of the electrons: f3 t h = 2 x 
1(T 3 T 4 1/2 . A similar expression is valid for losses due to the 
ionisation of neutral material. Using a composition appro priate 
to the interstellar medium, Mannheim & Schlickeiser (1994) 
find 
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where /3oc is the speed of an electron in the ground state of the 
hydrogen atom: /3o ~ 0.01. 




(7-l)Mc 2 (eV) 



frames (Kirk & Schneider |l98S|; Malkov & Volk |l995| ) and the 
self-excited waves are not damped by thermal ions, we can 
expect that k\ > 1. In this case, the constraint is not restrictive 
for the interstellar medium. 



10. Conclusions and implications for observations 

With respect to observational tests of shock acceleration in su- 
pernova remnants, the limit imposed by upstream wave damp- 
ing on the maximum particle energies is significant, but not 
serious, for the observability of supernova remnants in gamma- 
rays using the atmospheric Cherenkov technique. The limit can 
of course be circumvented by locating the particle acceleration 
and the gamma-ray production target in different regions or 
phases. In a clumpy medium, one could imagine accelerating 
in the low-density interclump phase with the dense clumps 
behind the shock providing the target material. Or one can 
consider the possibility of a SNR exploding near a molecular 
cloud (Aharonian et al. 1994 ). In addition one should consider 
the pre-ionisation of the upstream medium by the soft X-ray 
and UV radiation from behind the shock, an effect which will 
also reduce the effect of ion-neutral damping. 

One other possibility for obtaining observational evidence 
for shock acceleration in SNRs is th e use of Balmer di agnos tics 



in non-radiative shocks (Raymond 1991 



Smith et al. 



1994) to 



probe the structure of the cosmic ray precursor. This technique 
essentially determines the upstream plasma temperature and 
velocity averaged over the ion-neutral charge exchange length 
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which is numerically 
U 
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Fig. 7. The loss rates for ionisation (tzli dashed line) and Coulomb 
(*q . dotted line) losses compared to the acceleration rate for dif- 
fusive shock acceleration (tacc : solid line). Curves are plotted for 
parameters appropriate to a supernova occurring in a dense stellar 
wind: n H = 10 9 cm" 3 , Xi = 0.99, B = lmG. The spatial diffu- 
sion coefficient upstream and downstream of the shock is assumed 
equal to the Bohm value. The shock speed is 10 4 kms — 1 and the 
compression ratio is taken to be 4. 



This has to be compared with the characteristic scale associ- 
ated with particles of energy E, kb(E)/U. The two are equal 
at 

(lTw) = °'° 4 (lO 3 kms- 1 ) (lcm-0 



\lfxG J 



(58) 



The Coulomb and ionisation losses are roughly constant 
for particles of speed below /3 t hC and /3qc respectively and then 
fall as /3~ 3 at higher velocities, see Fig. ^. Comparing this be- 
haviour with that of t~^ c it is clear that the ionisation and 
Coulomb losses will not suppress the shock acceleration at any 
energy if the acceleration exceeds the loss rate for particles of 
speed /3 t h c and j3oc respectively. This is the case if 
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(55) 



Provided particles are injected at a speed several times the 
ion thermal velocity, so that the distribution function can be 
approximately isotropic in both the upstream and downstream 



The balmer diagnostic technique can in principle detect that 
part of the precursor structure produced by particles above 
this energy. 

In conclusion, the main result of this paper is that a more 
detailed analysis of the constraints on particle acceleration in 
dense media does not s uppor t the pessimistic view expressed 
by Draine and McKee ( |l993j ) "shocks that are optically ob- 
servable, either as nonradiative shocks in a partially neutral 
medium or as radiative shocks, are generally unable to accel- 
erate particles to extremely relativistic energies" . As far as the 
constraints imposed by upstream ion-neutral wave damping, 
and ionisation and Coulomb losses are concerned optically vis- 
ible SNR shocks should be capable of accelerating cosmic rays 
either to the revised cutoff ([Iq), or, if the condition (^) is sat- 
isfied to the maximum energy allowed by the shock geometry 
and age. 
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